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Abstract 

We present numerical and experimental results for the development of islands of sta- 
bility in atom-optics billiards with soft walls. As the walls are soften, stable regions 
appear near singular periodic trajectories in converging (focusing) and dispersing 
billiards, and are surrounded by areas of "stickiness" in phase-space. The size of 
these islands depends on the softness of the potential in a very sensitive way. 

Key words: billiards, soft-wall, atom-optics, cold atoms 
PACS: 05.45.Ac, 32.80.Pj, 42.50.Vk 



1 Introduction 



Kolmogorov-Arnold-Moser (KAM) theory provides a framework for the un- 
derstanding of an integrable Hamiltonian system which is subject to a per- 
turbation [1]. As the strength of the perturbation is increased areas of phase- 
space became chaotic and eventually only islands of stability, surrounded by a 
chaotic "sea", survive. The opposite question is of interest as well: What will 
happen to a completely chaotic system when a perturbation is applied? will 
it remain chaotic or will islands of stability appear? The appearance of KAM 
islands is of great importance: First, for initial conditions falling inside the 
island, there will be finite range oscillations in the the tail of the correlation 
function. In addition, these islands are usually accompanied by areas of " stick- 
iness" where particles belonging to the chaotic region of phase space spend 
an anomalous amount of time, and hence they affect the temporal correlation 
function even for initial conditions in the chaotic " sea" . Indeed, it was argued 
that this stickiness can cause a power-law decay of correlations, instead of an 
exponential one, and thus affect the transport properties of the whole system 
even for a relatively small island size [2]. 
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For the billiard system (i.e. a particle moving freely in a bounded region and 
reflecting elastically from the boundary), one such perturbation is making 
the walls "soft", as opposed to the infinitely steep potential of the ideal case. 
Since the billiard is a widely used paradigm for understanding the microscopic 
foundations of statistical mechanics, and noting that real classical molecules 
move in smooth potentials, we conclude that the study of soft-wall billiards 
is an important step towards sustaining these foundations and, in particular, 
Boltzmann's ergodic hypothesis. 

For a certain kind of dispersing billiards it was theoretically proven [3, 4] that 
when the wall becomes soft, an island appears near a singular, tangent tra- 
jectory (see details in section 3). More recently, and partly inspired by the 
numerical results presented herein, the appearance of an additional island, 
this time near a corner, was also shown[5]. Experimentally, some conjectures 
were made about the role played by the softness of walls in semi-conductor 
quantum dots. It was considered to be responsible for the algebraic tail in the 
distribution of the length of trajectories, which gives rise to the fractal nature 
of conductance fluctuations[6,7,8,9]. Observations of the formation of stable 
islands in a stadium shaped soft atom-optics billiard were reported in Ref. 
[10]. 

In this article, we present additional results of numerical and experimental 
studies of soft wall atom-optics billiards. Three unifying concepts are found in 
our results: First, islands of stability arise around singular trajectories (tan- 
gencies, corners) as predicted by Refs. [3,4,5]. Next, the size of these islands 
depend on the wall's softness in a very sensitive way. Finally, "stickiness" 
regions surround the stable islands and affect the decay properties of the bil- 
liards. In section 2, we present our numerical and experimental studies of 
the stability induced in an atom-optics billiard of the Bunimovich typefll] 
by making the walls soft. In section 3, we present results of numerical sim- 
ulations for the dynamics of atoms inside totally dispersing billiards, which 
exhibit softness induced stability around singular trajectories. 



2 Islands of stability in converging billiards 

The atom-optics billiard is a recently developed experimental system in which 
billiard dynamics can be studied[12,13]. As described in detail in Ref. [13], 
the billiard is realized by the use of a laser beam which is rapidly (100 kHz) 
scanned using two perpendicular acousto-optic scanners (AOSs). The laser 
frequency is tuned above the atomic resonance (the D2 line of 85 Rb in our 
experiment), hence it applies a repulsive force on the atoms. By controlling 
the deflection angle of both AOSs, we create the required billiard shapes which 
confine the atoms in the transverse direction. The instantaneous potential is 
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Fig. 1. Numerical simulations for the decay, through a small hole, from a tilted 
Bunimovich stadium (solid lines) and a circular billiard (dotted lines) with the 
parameters specified in the text, for several values of the softness parameter. Both 
billiards have the same area and hole size. The numbers adjacent to the lines indicate 
the value (in microns) of the softness parameter wq . The simulations for the stadium 
show a monotonic slowing down in the decay and increase in stability as the walls 
are made softer. For the circular billiard, there is almost no change in the decay 
curve as a consequence of the higher softness. The dashed line shows exp(-t/r c ), 
where r c is the escape time calculated for the experimental parameters. 

given by the dipole potential of the laser gaussian beam: 



where the curve (xo (t) ,yo (£)) is the shape of the ideal billiard, along which 
the center of the Gaussian laser beam scans, wo is the laser beam waist and 
Uq is the instantaneous potential height. The time averaged potential is kept 
constant for different values of Wq by changing the laser power. In such way, 
Wq serves as the softness control parameter, and is experimentally controlled 
by the use of a a telescope with a variable magnification, located prior to the 
AOS's such that w can be changed without affecting the billiard's size and 
shape. Fast enough scanning of the beam results in an effective time-averaged 
potential wall [14] . Photon scattering from the laser beam is greatly suppressed 
due to the large detuning, and collisions between the atoms are very rare for 
the range of densities used, hence motion of the atoms between reflections 
from the light walls can be regarded as strictly ballistic [1 3] . 



U(x,y,t) = U exp [-2 ((. 



x — 



x (t)) 2 + (y-y (t)) 2 ) /w 2 ], 
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Fig. 2. Poincare surface of section for monoenergetic atoms confined in a soft 
(wq = 24^m) tilted-stadium shaped atom-optics billiard with parameters as in Fig. 
1. The Figure shows ~ 9% of the total phase space area. A phase space area corre- 
sponding to chaotic trajectories is marked "A", one of period 2 trajectories inside 
an island "B", an area of higher periodicity trajectories "C" and one of "sticky" 
trajectories "D". 

The loading scheme of cold atoms into the billiard was described in [10]. 
Briefly, laser cooled 85 Rb atoms are loaded from a magneto optical trap into 
a red-detuned one-dimensional optical lattice which is laterally shifted (250 
/iin) from the billiard's location. The atoms are then transferred from the lat- 
tice into the billiard by pushing them with a pulse of a strong on- resonance 
beam which is perpendicular to both the lattice and the billiard beams. Simul- 
taneously with the pushing, the lattice beams are adiabatically switched-off, 
resulting in additional cooling[15]. Further reduction of the radial velocity 
spread and especially a decrease in the number of very slow atoms is achieved 
by capturing only a central velocity group of atoms into the billiard, which 
is turned on at the proper time after the push. Typically, 3 • 10 5 cold 85 Rb 
atoms are loaded into the billiard, with a typical velocity of ~ 100 mm/s and 
a velocity spread of ~ 20 mm/s. The cooling in the longitudinal direction 
ensures that the system can be approximated as a two dimensional one. 

As a way to probe the nature of the dynamics of the atoms confined by such an 
"atom-optics" billiard we measure the decay of the number of confined atoms 
through a small hole on the boundary. A purely exponential decay is expected 
for chaotic motion[16, 17,18], with a decay time given by r c = nA/vL, where 
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Fig. 3. Typical trajectories for atoms confined in a soft (wq = 24/tm) tilted-stadium 
shaped atom-optics billiard with parameters as in Fig. 1. (A) A chaotic trajectory, 
corresponding to the area marked "A" in Fig. 2. (B) A periodic trajectory of period 
2 and belonging to the central island, marked "B" in Fig. 2. (C) A high periodicity 
trajectory (Smaller islands surrounding the central one and denoted "C" in Fig. 
2). (D) A "Sticky" trajectory, spending in the island's vicinity a relative long time 
("D" in Fig. 2). 

A is the billiard's area, v the atomic velocity and L the length of the hole [16]. 
A slower algebraic decay usually characterizes a stable or partly stable phase 
space. In our system, the hole is produced by switching off one of the AOSs for 
~ 1 /is every scan cycle, synchronously with the scan. The number of atoms 
remaining in the trap is measured using fluorescence detection [14]. The ratio 
of the number of trapped atoms with and without the hole, as a function of 
time, is the main data of our experiments. 

First, we study a tilted Bunimovich stadium (see Fig. 3), which is chaotic and 
ergodic for the ideal hard- wall case. We used a "tilted" stadium, to avoid the 
long segments of regular motion associated with "bouncing ball" trajectories in 
the original stadium[19]. It is composed of two semicircles of different radii (64 
/mi and 31 /mi), connected by two non-parallel straight lines (192 /im long). 
The results of numerical simulations for the decay in the number of confined 
atoms, through a small hole in the boundary, are presented in Fig. 1. For the 
stadium with w = 12 /im a nearly exponential decay is seen, as expected 
for the ideal hard-wall billiard. As the walls are made softer, a monotonic 
slowing down in the decay is observed. We compare this behavior with that 
of a circle with the same area and hole size, which exhibits integrable motion. 
In the ideal hard-wall case the existence of nearly periodic trajectories results 
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Fig. 4. Experimental results for the decay of cold atoms from a tilted-stadium 
shaped atom-optics billiard, with two different values for the softness parameter: 
wq = 14.5/xm (o), and wq = 24fim (+). The hole is located inside the big semicircle. 
The smoothening of the potential wall causes a growth in stability, and a slowing 
down in the decay curve. 

in many time scales for the decay through a small hole on the boundary, and 
then yield an algebraic decay [16]. Simulations for the decay from a circular 
billiard (dotted lines in Fig. 1) show that the decay is almost unaffected by 
the change in softness parameter, as expected. 

Figure 2 shows the Poincare surface of section obtained from the results of 
numerical simulations for classical trajectories of Rb atoms inside the tilted- 
stadium billiard, with wq = 24/im. For clarity, we assume a monoenergetic 
ensemble (with v = 120 mm/s), a two-dimensional system, and no gravity. 
However, including the experimental velocity spread, three dimensional mo- 
tion and gravity have only a small effect on the results. In analogy to our 
experimental setup, we use for our simulation a scanning gaussian potential, 
with a scan frequency of 100 kHz [20]. The dimensions of the billiard are equal 
to the experimental ones. The Poincare surface of section shows v x versus x 
at every trajectory intersection with the billiard's symmetry axis (y = 0), 
provided that v y > 0. The smoothness of the wall results in a stability re- 
gion that appears around the singular trajectory which connects the points 
where the big semicircle joins the straight lines. This fact is also confirmed 
by analyzing some typical trajectories (see Fig. 3). Additional islands (smaller 
than the central one) appear around it, and correspond to trajectories with 
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Fig. 5. Experimental results for the decay of cold atoms from a circle shaped 
atom-optics billiard, with different values for the softness parameter: wq = 14.5/xm 
(o), w = 18.5/xm (+), wo = 21/xm (□), w = 23.5/Ltm (x). 

higher periodicity, like the one shown in Fig. 3C. Around these islands there 
is a large area of "stickiness", where the trajectories spend a long time. Such 
a "sticky" trajectory is presented in Fig. 3D. The exact structure of the island 
and its vicinity depends on the softness parameter wq in a sensitive way [10]. 
In general, the size of the island and the stickiness region increase with the 
softness parameter. The slowing down in the decay is attributed mostly to the 
sticky trajectories, which in general cover a phase space volume greater than 
that of the island itself (see Fig. 2 of [10]). 



In Fig. 4, experimental results for the decay from a tilted stadium shaped 
atom-optics billiard are presented, for two different values of the softness pa- 
rameter (w = 14.5 jum and w = 24 /zm). The experimental results confirm 
that the soft wall causes an increased stability, and a slowing down in the de- 
cay curve, in good agreement with the numerical results. We also observed[10] 
that when the hole includes the singular point where the semicircle meets the 
straight line, no effect for the change in wq is seen. In this case the stable island 
that is formed around the singular trajectory and the region of "stickiness" 
around it (see Fig. 3) are destroyed by the hole. By moving the hole we can 
thus experimentally map the location of the stable island. 

Similar decay measurements for a circular atom-optics billiard (see Fig. 5) 
showed no dependence on w in the range 14.5 — 24 fj,m. We were not able to 
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Fig. 6. Poincare surface of section calculated for monoenergetic atoms confined in the 
soft dispersing billiard described in the text with £ = 40 and wq = 24 /jm. The 
Figure shows ~ 17% of the total phase space area. Phase space area corresponding 
to chaotic trajectories are marked "A" , those of period 2 trajectories inside an island 
"B", areas of higher periodicity trajectories "C" and "sticky" trajectories "D". 

see any dependence on the hole position either. 



3 Islands of stability in dispersive billiards 



We studied the effect of a soft potential wall also in a dispersing billiard. 
The billiard that we studied was formed from the intersection of four arcs 
with different radii and it had a tangency, namely, a specific periodic orbit is 
tangent to the billiard's boundary at one point (see Fig. 7). In this billiard, 
the periodic orbit had a period of two. More generally, we consider a family 
of such billiards, with a variable parameter (, which is the vertical distance 
between the lower arc and the two-periodic orbit. For ( = the periodic orbit 
is tangent to the lower arc. For ( > the lower arc is below the periodic orbit, 
while for ( < it is above and actually intercepts the orbit. A billiard with a 
two-periodic tangent orbit was chosen, to increase the size of the island. Every 
hard-wall billiard of this family is fully chaotic. This specific billiard shape was 
motivated by a theoretical study which predicts the formation of an elliptic 
island around such a tangency [4]. 



Fig. 7. Typical trajectories in a soft dispersive billiard, with the parameters as in 
Fig. 6, and softness parameter wo = 24 fim. (A) A chaotic trajectory. (B) A periodic 
trajectory of period 2. (C) A high periodicity trajectory. (D) A "Sticky" trajectory. 

In Fig. 6, the Poincare surface of section in the vicinity of the two-periodic orbit 
is presented, for a billiard with £ = 40 /im, and Wq = 24 /im. An elliptic island 
is clearly observed. Typical trajectories for this case are shown in Fig. 7. In Fig. 
7A, a typical chaotic trajectory is shown, corresponding to the areas in phase 
space marked as "A" in Fig. 6. A periodic trajectory of period 2 (corresponding 
to the central island in Fig 6, marked as "B") is shown in Fig. 7B. Also shown 
in Fig. 7 are a high-period periodic trajectory and a "sticky" one. We repeated 
these calculations for billiards with different values of ( and wq, and plotted 
the size of the island as a function of these two parameters in Fig. 8. It can be 
seen that the island becomes observable above some softness of the potential, 
and that its size initially increases with an increased softness. However, further 
softening reduces the size of the island. In general, higher values of ( require 
higher values of wo to induce the stability. These results, although obtained 
for a strongly perturbed case, are in agreement with Refs. [3,4]. Some intuition 
can be gained by noting that the stable trajectory is slightly reflected by the 
lower arc. Since the stable trajectory is almost tangent to the lower arc, its 
perpendicular velocity is very low, and a small potential is already sufficient to 
reflect it. This small potential is provided at a certain distance from the lower 
arc, depending on wo. If ( is increased, the trajectory will "feel" the lower arc 
(and hence will be stabilized) only for a larger value of w . However, if w is 
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Fig. 8. Island size in the soft wall dispersing billiard, as a function of the softness 
parameter, wo, and the shape parameter, £. (White - no island with an area greater 
than 1.6 x 1CT 3 of the total phase space area, black - a large island, > 4 x 1CT 2 of 
the total phase space). 

further increased, the potential at the trajectory would become too high, thus 
"blocking" it and destroying the stable island. In additional simulations we 
observed that the center of the island moves towards larger y values when wq 
is increased, in good agreement with the above intuitive model. 

We observed also another kind of stable islands in these billiards, around a 
different type of singular trajectories, i.e. trajectories that are supported by the 
corners of the billiard. An example for such a trajectory, and a sticky trajectory 
in the area surrounding this new island, are shown in Fig. 9. Our observations 
are consistent with the theoretical proof for the existence of this large islands, 
provided in Ref. [5]. There it was shown that, for a small perturbation, its 
existence depends both on the geometry of the corner and on the smooth 
potential behavior at the corners. Our results show the existence of this island 
also for a strong perturbation. In all cases where an island was found, we also 
found a sticky region near it with atoms spending arbitrary long times near 
the island. 

We performed decay simulations for the scattering billiard, by opening a hole 
on the boundary, and counting the number of remaining atoms as a function 
of time. Here, as in the case of the Stadium billiard, an exponential decay over 
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Fig. 9. Corner trajectories. (A) A periodic trajectory of period 2. (B) A "Sticky" 
trajectory. 

3 orders of magnitude was observed for wq = 8 jum. In cases where an island 
develops, a slowing in the decay was observed when the hole did not include 
a reflection point of the periodic orbit. This slowing down in the decay was 
similar to that for the stadium billiard of Fig. 1. 



4 Summary 

In this paper we showed, numerically and experimentally, that stable islands 
are created in the phase space of chaotic billiards, either dispersive or focus- 
ing, when the walls are softened. Using cold atoms in atom-optics billiards, an 
island of stability was demonstrated in a tilted Bunimovich stadium billiard, 
around the periodic trajectory that connects the points where the big semicir- 
cle joins the straight lines. For a dispersive billiard, we confirmed numerically 
the existence of a stable island around a periodic orbit which is tangent to the 
billiard's wall. In addition, we found numerical evidence, for the first time, for 
the formation of a stable island near a "corner" trajectory. In all three cases, 
the elliptic islands appear around trajectories which are singular in the "ideal" 
infinitely hard wall billiard. The orbits are singular in the sense that they do 
not have a smooth limit when approaching the singularity from both sides. 

Our numerical results show that the appearance of a stable island in a soft- 
wall billiard is in general accompanied by areas of "stickiness" surrounding it. 
Although the size of the stable island is a sensitive function of softness, not 
always showing a monotonic behavior, we confirmed numerically and experi- 
mentally that the "sticky" trajectories modify this behavior and, in general, 
introduce a monotonic slowing down in the decay of the number of trapped 
particles, as the walls are made softer. 

An interesting question in this context is wether the existence of a singularity 
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in the hard-wall billiard is a necessary and sufficient condition for stability in 
the soft-wall one. More generally, one can ask to what extent is the relation 
between the existence of a singularity in a chaotic system and the emergence of 
stability when this system is perturbed, a universal feature. From the practical 
point of view, the softness of the wall can be used as a control parameter for the 
structure of phase space, in particular in the context of controlling chaos[21]. 
Finally, we note that since stable and "sticky" regions follow the softness of 
the walls, and since any physically realizable potential is inherently soft, the 
physical realization of a truly chaotic and ergodic billiard is questioned. 
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